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We present a covariant decomposition of Einstein's Field Equations which is particularly suitable 
for perturbations of spherically symmetric - and general locally rotationally symmetric - spacetimes. 
Based upon the utility of the 1+3 covariant approach to perturbation theory in cosmology, the semi- 
tetrad, 1+1+2 approach presented here should be useful for analysing perturbations of a variety of 
systems in a covariant and gauge-invariant manner. Such applications range from stellar objects 
to cosmological models such as the spherically symmetric Lemaitre-Tolman-Bondi solutions or the 
class of locally rotationally symmetric Bianchi models. 
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, Tetrad formalisms in general relativity have played a pivotal role in its development as well as our understanding 
of the subject. These range from the complex null tetrad of Newmann and Penrose, to the 1+3 approach of Ehlers, 
Ellis and others, which includes both a full tetrad approach as well as a partial 'covariant' approach where only one 
timelike tetrad vector is chosen (see [l| for a review and references). These techniques formulate the equations of 
i— 1 ' general relativity as first-order differential equations in the physical curvature and dynamic variables of the covariant 
derivatives of the tetrad vectors, as opposed to the more usual coordinate approach involving second-order partial 

' differential equations in functions appearing in the metric. The differential operators which appear are convective 
derivatives along the tetrad vectors as opposed partial derivatives with respect to particular coordinates. Much of their 

1 utility arises in spacetimes with special symmetry. For example, the 1+3 covariant approach is perfect for cosmology 
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I. INTRODUCTION 



because it covariantly factorises out the essential coordinate - time - leaving all the background field equations as 
— i • covariant scalar equations. Under perturbations all 3-vectors and tensors (which must vanish in the background due 
to homogeneity and isotropy) become gauge- invariant first-order quantities making a Fourier analysis easy [l], Q • 

We formulate here an approach which involves a semi-tetrad: we keep the timelike threading vector field of the 1+3 
approach and introduce one spatial vector. The remaining two dimensions are left untouched, rather like the '3' in 
the 1+3 approach. Indeed the formalism presented here may be considered as half-way between the 1+3 tetrad and 
covariant approaches. A similar approach has been discussed before in [1,0, Hill, an d we expand on this considerably 
here by presenting the full system of 1+1+2 equations. 

It is expected that this approach may find use in perturbations of spacetimes with a preferred spatial direction 
at each point - so-called locally rotationally symmetric spacetimes Q- These include the spherically symmetric 
Lemaitre-Tolman-Bondi models, many classes of Bianchi models, as well as forming the background for most stellar 
models. In this paper we provide the algorithm of how to calculate gravitational perturbations in any LRS spacetime 
in a covariant and gauge-invariant (GI) way. 

As an example of it's utility, such a covariant perturbative scheme was applied to the Schwarzschild solution in [8[. 
?3 i Despite being a well understood problem, it was shown using the 1+1+2 approach how both the axial and polar 
degrees of freedom may be unified into a single transverse traceless tensor which obeys the tensorial form of the 
Regge- Wheeler equation [1, Q 

-W ab + & ab + AW ab -tfW ab + 6 2 W ab = 0, (1) 

where the Regge- Wheeler tensor W ab is a gauge- and frame-invariant TT tensor, defined in ]|] (other variables are 
defined below), and and 6 are time, radial and angular derivatives, respectively. This tensor contains in compact 
form the curved space generalisation of the two flat space GW polarizations h+ and h x (see also [ll[ for an 
extension of this work) . The approach here also been used to studyscalar and electromagnetic perturbations of LRS 
spacetimes, and generalised Regge- Wheeler equations were found [I3,ll3j|. Furthermore, it has been used to study the 
interaction of magnetic fields and gravitational waves around a black hole - a process which produces electromagnetic 
radiation mirroring the gravitational waves [To| . 
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In Sec. [IT] we discuss the 1+1+2 approach in full generality, and then in Sec. Ill we discuss the perturbation 
procedure for LRS spacetimes, before sumarising in Sec. IV. 



II. FORMALISM 



In the 1+3 approach, a timelike threading vector field u a (u a u a = —1) is introduced, representing the observers' 
congruence. Given this vector field, the projection tensor h b = g b + u a u b is introduced, which projects all vectors 
and tensors orthogonal to u a . Using h a b, any 4-vector may be split into a (1+3 scalar) part parallel to u a and a 
(3- vector) part orthogonal to u a . Any second rank tensor may be covariantly and irreducibly split into scalar, vector 
andprojected, symmetric, trace-free (PSTF) 3-tensor parts, which requires the alternating tensor e a bc = u d r]dabc Q]- 
Tensors of higher rank may be similarly split, but are rarely used (an important exception being cosmic microwave 
background physics [HHHI]). These are the fundamental quantities describing the spacetime, after the introduction 
ofu a . 

We now introduce another vector field and perform another split, but this time of the 1+3 equations. The '1+1 +2' 
decomposition we develop here has been partially studied before, mostly in the context of symmetries of solutions of 
the EFE 0,(3. It was introduced by [|| and further developed in pi la. [Til Il2l [l3| . However, there are importances 
differences with the work presented here. In the following we assume the 1+3 covariant split of the equations (as 
given in [l[, for example), with all tensors split into scalars, vectors and PSTF tensors with respect to u a . 

Take a unit vector n a orthogonal to u a : n a n a — 1, u a n a — 0, and define the projection tensor 

N b = h b -n a n b =g b + u a u b -n a n b , (2) 

which projects vectors orthogonal to n a (and u a ): n a N a b = = u a N a b, onto 2-surfaces (N a a = 2) which we refer to 
as the sheet. This is also the screen space of the null vector k a cx u a + n a . 

Any 3-vector ip a can now be irreducibly split into a scalar, 'J, which is the part of the vector parallel to n a , and a 
vector, lying in the sheet orthogonal to n a ; 

i; a = ^n a + y a , where # = ij} a n a , and * a ee N ab ^ b = *p a , (3) 

where we use a bar over an index to denote projection with N a b on that index. Similarly, any PSTF tensor, ipab, can 
now be split into scalar, vector and tensor (which are PSTF with respect to n a ) parts: 

V'afc = ip(ab) = * (n a n b - ±N ab ) + 2* (a n h) + * afc , (4) 

where 

* ee n a n b i> ah = -N ab i) ab , 
* a ee N a b n c ili be = * fi , 

* Qb ee ^ {ab} ee [N (a c N b) d - iiV afc A crf ) ^ cd . (5) 

We use curly brackets to denote the PSTF with respect to n a part of a tensor. Note that for 2nd-rank tensors in the 
1+1+2 formalism 'PSTF' is precisely equivalent to 'transverse-traceless'. Note also that h^ ab j — 0, = —n^ a n b ) = 

N a b — f^ah- 
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We also define the alternating Levi-Civita 2-tensor 

£ab = £ab c n c = u d i] dabc n c , (6) 

so that e a bn b — = e( M, and 

£ a bc = n a e bc - n b e ac + n c e ab , (7) 
s ab s cd = N a c N b d ~N a d N b c , (8) 
e a c e bc = N ab , e ab e ab = 2. (9) 

Note that for a 2-vector e a b may be used to form a vector orthogonal to ^ a but of the same length. 

With these definitions we may split any object into scalars, 2-vectors in the sheet, and transverse-traceless 2- 
tensors, also defined in the sheet. These three types of objects are the only objects which appear, after a complete 
decomposition. Hereafter, we will assume such a split has been made, and 'vector' will generally refer to a vector 
projected orthogonal to u a and n a , and 'tensor' will generally mean transverse-traceless tensor, defined by Eq. (J5|). 
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There are two new derivatives of interest now, which n a defines, for any object ip--- " : 

^ a ... b c - d = n e D e ^ a ... b c- d , (10) 
5^ a ... b c - d = NjNj ■■■N b 3 N h c ---N i d T>^ f ... g h - i . (11) 

The hat-derivative is the derivative along the vector field n a in the surfaces orthogonal to u a . This definition represents 
a conceptual divergence from 1+3 tetrad approach, in which the basis vectors appear on an equal footing [i.e., with 
V a rather than D a in Eq. (|10|)]. As a result, the congruence u a retains the primary importance it has in the 1+3 
covariant approach. (We choose to think of A = u a n b V a Ub = —u a u b V ' a nb as the radial component of the acceleration 
of u a , rather than the time component of h a .) The (^-derivative, defined by Eq. (fTl"]) is a projected derivative on the 
sheet, with projection on every free index. 

These derivatives then affect our projection tensor N a b and Levi-Civita tensor as follows: 

Nab = 2u( a u b )-2n( a h b )=2u( a Ab)-2n( a a b ), (12) 

N ab = -2n.( a n b ), (13) 

ScNab = 0, (14) 

e a b = -2u [a s b ]cA c + 2n [a e b ] c a c , (15) 

4b = 2n [a e b ] c a c , (16) 

S c e ab = 0. (17) 

We now decompose the covariant derivative of n a orthogonal to u a into its irreducible form: 

B a n b = n a a b + \<j>N a b + £,£ab + Cab, (18) 

where 

a a = n cr D c n a = h a , (19) 

<P = S a n a , (20) 

£ = y ab 6 a n b , (21) 

Cab = 6{ a n b }- (22) 

We may interpret these as follows: travelling along n a , (j> represents the sheet expansion, Cab is the shear of n a 
(distortion of the sheet), and a° its acceleration, while £ represents a 'twisting' of the sheet - the rotation of n a 0]. 
The other derivative of n a is its change along u a , 

n a = Au a + a a where a a = h a and A — n a u a - (23) 

The new variables a ai (f>, £, Cab an d cn a are fundamental objects in the spacetime, and their dynamics gives us 
information about the spacetime geometry. They are treated on the same footing as the kinematical variables of u a 
in the 1+3 approach (which also appear here). 

For any vector fy a orthogonal to n a and u a (i.e., ^ a = 5 ,a ), we may decompose the different parts of its spatial 
derivative: 

B a ^b = -n a n b y c a c + n a ^ b - n b [±0tf a + (£e ac + £oc) * c ] + S a 9 b . (24) 
Similarly, for a tensor ^ a b- ^ab = ^{ a b}, we have 

D a * bc = -2n Q n (b * c)d a d + n a * bc -2n (b \<j>^ c)a + ^ c) d {& ad + Cad) + S a ^ bc . (25) 

Note that for a scalar, we have D a ^> = ^>n a + S a ^>. 

We take n a to be arbitrary at this point, and then split the usual 1+3 kinematical and Weyl quantities into the 
irreducible set {9, A, ft, S, £, H, A a , S a , £ a , H a , £ a6 , £ ab , U ab ) using © and Q as follows: 

ii a = An a +A a , (26) 

uj a = ftn a + fi a , (27) 

a a b = S (n a n b - \N a b) + 2E( a n h) + E a6 , (28) 

Eab = £ {n a n b - \N a b) + 2£( a n b ) + £ab, (29) 

H ab = H (n a n b - ±Nab) + 2H (a n b) + H ab - (30) 
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The shear scalar, a, for example, may be expressed in the form 

a 2 = \o ab a ab = f E 2 + E Q E a + lE ab E ab . (31) 

Similarly we may split the fluid variables q a and 7r ab , 

q a = Qn a + Q a , (32) 
ir ab = II (n a n b ~ ^N ab ) + 2n (a n o) + Tl ab . (33) 

Having described the splitting of the 1+3 variables to obtain their 1+1+2 parts, and the introduction of the new 
1+1+2 variables corresponding to the irreducible parts of V a ro b , it only remains to apply this decomposition procedure 
to the 1+3 equations themselves, as well as the Ricci identities for n a . We give these equations in section HTT1 



A. Commutation relations 



In general the three derivatives we now have defined, ' ' ', '"' and 'S a ' do not commute. Instead, when acting on a 
scalar ip, they satisfy: 



8 a ij - N b (5 b ipy = -A a i> + (a a + E a - e ab Q b ) $ + (§0 - ±£) S a1 p + (E ab + fte ab ) 
8 a ^-N b {5^) = (^ a ~£ ab Sl b )i, + a a 4> + \(t>8 a ^ + {t ab + te ab )8 b 4>, 
5 a 5 b ip - 5 b 5 a il> = 2e ab (q4> - + 2a [a 5 b] ip. 

The commutation relations for 2-vectors ip a are 



(34) 
(35) 
(36) 
(37) 



-Aip a + (§0 + E) fa + (E 6 + e bc n c - eft) <5>a + A a (E 6 + £ 6c ^ c ) 4 + We a6 </> 6 , (38) 
- A^fc + (a„ + S a - e ac n c ) 4>- b + (|(9 - |E) (c5 Q ^ b + VaA) + H a €u^ c 
+ (E ac + tte ac ) {S c i, b + i/; c A) + 5 (4> a Qb - N ab ^Qc) 

- (^<f>N ac + £e ac + ( ac ) ^ c a b , (39) 
(E a - s ac n c ) ip b + a a ijj b + i0 (S a ip b - ip a a b ) + (C ac + £e ac ) (8 c ip b - ip c a b ) 

+7v Qb ^ c (in c + £ c ) - ip a (±n b + s b ) , (40) 

2e ao (^4 - ftfe) + 2 [(£0 - iE) 2 - l0 2 + |n + £ - i (/i + A)] ^ [o iV b]c 

-2V» [a [- {\0 - |S) (S 6 ] C + ^fc]c) + 3^ (Cb]c + fo]c) + 2 n fc]c + fyc] 

+27V [ac [- (§0 - IE) (E b]d + fte b]rf ) + 10 (C 6]d + £e b]rf ) + \K b]d + £ b]d ] V d 

+2 [- (E [QC + Qe [ac ) (E b]d + Qe b]d ) + (C [ac + £e [oc ) (C b]ci + ^6]d)] V> d - (41) 

These relations are more complicated for tensors. These last two equations in the case of scalars are the decomposition 
of the 1+3 commutation relation 



S a ip b -N a c N b d (S c M 



Sj b - N a c N^(S c ^ d ) 

SaSblpc - SbSa^Pc 



curl D a ip = 2ipu) a . 



(42) 



From Eq. ()37|) . we see that our sheet will be a genuine 2-surface in the spacetime (and, in particular, that the 
derivative S a will be a true covariant derivative on this surface) if and only if £ = ft = a a = 0. (Recall that the 1+3 
spatial metric h ab corresponds to a genuine 3-surface when uo a = 0.) Otherwise, the sheet is really just a collection 
of tangent planes. In addition, the two vectors u a and n a are 2-surface forming if and only if the commutator [u, n] 
in has no component in the sheet: that is, when Greenberg's vector 



E + e il b — ot 



(43) 



vanishes 

" see Eq. QgJ). 
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III. THE EQUATIONS 

Once the vector n a has been introduced it is possible, and necessary, to augment the 1+3 equations with the Ricci 
identities for n a ; without these we do not have enough equations to determine the new 1+1+2 variables. The Ricci 
identities for n a are 

Rabc = 2V [a V b] n c - Rabc d n d = 0, (44) 

where R a bcd is the Riemann curvature tensor. This third-rank tensor may be covariantly split using the two vector 
fields u a and n a , and gives dynamical equations for the covariant parts of the derivative of n a (namely a a , a a , <fi, £ 
and Cab) in the form of evolution equations, involving dot-derivatives of these variables, and propagation equations, 
involving hat-derivatives. In order to facilitate the calculation of these Ricci identities, which appear in the following 
section, we give here the expression for the full covariant derivative of n a in terms of the relevant 1+1+2 variables: 

V a n fc = ~Au a u b - u a a b + (E + ±0) n a u b + (E a - e ac Q c ) u b + n a a b + \^N ah + £s ab + Cab, (45) 

which may be inserted into Eq. (|44[l . The full decomposition of the covariant derivative of u a is 

V a u b = -u a {An b + A b ) + n a n b (|0 + E) + n a (E 6 + e bc fl c ) 

+ (E„ - e ac n c ) n b + N ab (±0 - ±E) + Sle ab + E ab , (46) 

which in turn implies the useful relation 

u a = (±0 + E)n Q + E Q + £(lb ft b . (47) 

We have now assembled all the tools necessary to provide the full system of equations for the 1+1+2 formalism. 
This consists of evolution equations, propagation equations, mixtures of both, and constraints. Formulae which are 
useful for splitting 1+3 equations are given in Appendix [Al 

A. Evolution equations 

We find evolution equations for the 1+1+2 variables 0,£ and Cab from the projection u a R abc . 
u a N bc R abc : 

4> = (§0 - E) (A - \<f) + 2£Q + 6 a a a + A a (a a - a a ) + (a a - A a ) (E a - e ab n b ) - ( ab Z ab + Q, (48) 

n .a~bc. td 

U 8 K a bc- 

i = (|E - i0) Z+(A-^)n+t(a a + A a ) [fi a + e ab (a b + E 6 )] + ±£ Qb 6 a a fc - ±£ CQ C b c E afc + (49) 

M C i? c {ab} : 

C{ab} = (|E - ±0) Cab + Q£c{a( b } c + (A - \<j>) E ab - ££ c{a E b} c - Cc{ Q E b} c + 5{ a a b } 

+ (.4 {a - a {a ) a b} - (A {a + a {a ) (E b} - £ 6}d n d ) - e c{a H b} c , (50) 

Then a 1+1+2 decomposition of the standard 1+3 evolution equations gives us the remaining evolution equations, 
which can't be found from R abc - 
Vorticity evolution equation: 

f2 = \e ab 5 a A b + Ai + tl{^- §0) +r> a (E a + a a ) (51) 

Shear evolution: 

S{ a b} = <5{a*4b} + A{ a A b } - E{ Q [E 6 } + 2a 6 }] - ^{a^f,} + -4Cab - (§0 + 3E) E ab — E c {- a E 6 } C - £ ab + in ab (52) 



B. Mixture of propagation and evolution: 



u a n b R abs = n a u b R abe : 

a- a -a- a = - {\<j> + A)a a - & ab a b + (\6 + S) (A a - a a ) + {\<f> - A) (E„ + e ab Q b ) 
-Z (s ab Z b - n„) + Cab (-a b + S b + e bc n c ) + \Q a - e ab H b , 

u a n b u c R abc = -n a u b u c R abc : 

A-\0-t = -A 2 + (§0 + S) 2 - 2a a S a + S a S a - ft a ft a - a a -4 a + e ab a a Q b 
+ i(^ + 3p-2A)+£-in 

Raychaudhuri equation: 

A - 8 = -S a A a -(A + 0)A+(a a - A a ) A a + \6 2 + §S 2 - 2Q 2 + 2S a S a 
-2fl a n a + S ao S ab + | ( M + 3p) - A 

Vorticity evolution: 

+ \e ab A b = - (f + |S) Q a + n (S a - a.) + |64a + ±e a6 [--4a b + S b A - ^A b ] 
-\e ab ( bc A c + S ab fi b 

Shear evolution: 

S - lA = \(2A-(j>)A-(l6+ |S) S — |fi 2 — i<5 a ^ a + S a [2a a - iS a ] 
-\A a [2a a - A a ] + i0 a a + iS ab S ab - 5 + |n 

S a - i Aa = \5 a A +(A- \<t>) A a - (f + |S) S a + i.4a a - §Sa a - Qft a 
-| (£e«6 + Cab) A b + S ab (a" - S b ) - £ a + in a 

Energy conservation: 

fi + Q = -s a Q a -6(fj, + p)-(<t> + 2A)Q- §sn + (o„ - 2A) a - 2s Q n a - s ab n a ' 

Momentum conservation: 

Q+p + fl = -S a U a -{\<t> + A)U- (§0 + S)Q- (M + f>M 

+ (a„ - S a + e ab f! b ) Q a + (2a a - X) n a + Ca h n ab 



Qa + fla = S a p + \5Jl - S b Il ab - Q (a„ + S a + e afc fi b ) - §ITa a - (§0 - iS) Q a + ne ab Q b 

- (§</> + a) n a + ee ao n b - + p - in) A a - s ab Q b - c ao n b + n afc (a b - A b ) 

Electric Weyl evolution: 

£ + ill + \Q = +e ab 5 a H c + \5 a Q a + (§E - 0) £ - \ (§0 + §£) n + § (\<j> -2A)Q + 3£W 
-i ( M + p) S + (2a a + S a - e ab tt b ) £ a + (a„ - iS a - i £afe Q 6 ) n a 
+i (o„ + A) a + 2e ab AW c - S ab (£ ab + in ab ) + e ab H bc C c 



£- a + \e ab H b + ili a + \Qa = \e ab 5 b H + \e bc 5 b U c a - \8 a Q -\{»+P-le+ in) S a 

+| (£ + in) e ab n b - \QA a + \neabA b -!{£ + |n) « a 

+i (i0 - ^) Q a + \ie ab Q b -\{\e + is) n a - ir>e ab n b 

+is ab (3f b - in 6 ) + i (3£ o6 - in Qb ) s b - (5 ah + in ab ) (a 6 + i e bc r> c ) 

+ 2Ca6 (e C W C + Q ) — Hate C ^4 C 
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£{ab] ~ £ C {cfib} + l^{ab} = ~£c{ a 5 C Kb} ~ 5<5{aQb} - \ + P + 3£ - 5II) E ab - ^QC,ab - f He c { a C b } C 

- (9 + §£) £ ab + fte c{a £ b} c - ($0 - |£) n ab + ifie c{a n b} c + £W a b 
+ (\<j> + 2A) e c{a H b] c - A {a Q b} - (a {a + ±£ c{a n c ) (2£ b} + H b} ) 
+E {a (3£ b} - |n b} ) + 2e c{a H b} (a c - A c ) 

+E c{a (i£ b} c - in b} c ) + e c{a H b}d C d (64) 

Magnetic Weyl evolution: 

H = -e ab 8 a £ b + y ab S a U b - 3££ + (9 + §E) W + SIQ + §£II - 2e ab A a £ b 

+ (2a a + E Q - e ab n b ) U a -\ (fi„ + e ah S b ) Q a - E Qb H ab - \e ab £ b % a c (65) 

W 5 — \e ab £ b + \e a bR h = ~j£ a bS b £ + ^e ab S b U — ^e bc S b £ c a + je bc S b H c a 

+ jHT. a + \Qe ab Y, b + lQfl a + \Ue ab ^ b - f £e ab A b - \Ua a 

+ § {£ - in) e ab a b - \i£ a + {\<f> + A) e ab £ b + (§£ - 0) W a - ifte Q bH b 

+ §ftQ a - |Se ab Q h + |en a - i0e ab n b + E ab (§W b + i £ bc Q c ) 

+ f £ Qb C 6c (£c - §n c + f A) + H ab {a b + §E b - ie bc f! c ) (66) 



^{ab} + £ C {a£ b } C - l £ c{ a n b } C — £c{a<5 C £h} - 2 £ e{a<5 C n b } - §7YE ab + i|(3£ c { Q E b } C 

+| (£ - in) e c{a C b} c - i£ ab - {\<j> + 2A) e c{a £ b} c 

- (9 + §£) "H ab - rie c { a 7Y b } C + ^n ab + i(/)e c { a n b j c 

+E{ a (3H b } - e b } c Q c ) + ^{a (|Qb} _ £b} c H c ) - 2a^ a H b ^ 

+£{ Q 2e b } c (a c + yl c ) - n{ a e b } c a c + 3E c { a H b } C - e c { a C cd (£ b }d ~ \^b}d) (67) 

C. Propagation equations 

Propagation and constraint equations are formed from either projecting i? abc as indicated, or from projections of 
the 1+3 constraint equations, denoted Cj, as given in 
n a N bc R abc : 

$ = -\4> 2 + 2C 2 + {\0 + E) (f A - S) + S a a a - a a a a - UC" + 2e ab a a n b - E a E a + n a il a 

-|( M + A)-in-f, (68) 

n a e bc i? abc : 

I = + (§0 + E) + ie Qb 5 a a b + ie Qb E a a b + (£o„ + a.) il a (69) 

C{ab} = -<ft(ab - C{ a Cb}c + <5{aOb} - a{ a Ob} + 2a{ a £b}c^ C - ^{ A} - E{ a E b } 

+ + E) E ab - in Qb - £ ab , (70) 

Shear divergence (Ci) a n a : 

t-\§ = -l^~2Zri-5 a Z a -e ab 8 a n b + 2Z a a a -2e ab A a n b + Z ab e b -Q (71) 
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and (Ci) 5 : 

±- a -e ab n b = ±5 a X + l8 a 8~e ab 5 b n~^ a + & ab Y, b -& a + (^ + 2A)e ab n b 

-§Sa + ne ab [a b - 2A b ] - S b Z ab - ( ab Y, b + S ab a b + e ab C bc r> c - Qa (72) 

Vorticity divergence equation (C2): 

n = -5 a n a + (A~ci>)n + (a a + Aa)n a (73) 

{C3){ab}- 

S{ ab } = <S{ a S b } - £ c { a 5 C VL b } ~ \4>Y, ab + ££ c { a £ 6 } C + §SCab - ^£ C {aCfc} C 

-2S{ a a h j - 2e c { a ^l c ri b } - E c { a C 6 j. c - e c { a Ti h y (74) 
Electric Weyl Divergence (C4) a n Q : 

f - |/t + |n = -(5 a £ a - i<5 a n a - 10 (f + |n) + (|s - |6>) q + 3fiH + (2£ a + n a ) a a 

+ iS a g Q + 3n a H a - f e Qb O a Q b + e ab X ac H b + (Sab + ^Tlab) C b (75) 

(Ci) a : 

£a + |ft a = |<5 a £ + |(5 a /x + i<5 a n - 5 b £ ab - i(5 b n ab + iQS a + He ab S b - f wn a - §Q£ Qb r* b 
-| (£ + |n) a Q - §0 (£ a + in a ) + & ab (£ b + in 6 ) + 3ftW Q - Se ab H b 

- (§0 + is) Q a + ^e ab Q b + is ab g b - Cab (£ b + ^n b ) + (Sab + ^n Qb ) a b + m a b^ b (76) 

Magnetic Weyl divergence {C^) a n a : 

n = -8 a n a - \e ab 5 a Q b - \<m - (3£ + a* + p - §n) n - oe 

+2H a a a - 3tt Q (£ a - in Q ) + C a6 H ab - e ab S a c (£ bc + in bc ) (77) 

(C 5 )a. 

Ua-\eabQ b = \SaU-5 b Ua b ~\e ab 5 b Q~l(£ + in)e ab S b - (-§£ + M + p + in) O a 
-|Wo + iQ£ ah a b - 3Q£ a + §Se ab £ b - \4>H a + C~e ab H b 

-\iQ a + \4>e a bQ b + |on a + |Ee ab n b 

+n ab a b - Q ab U b - 3 (Sab - in ab ) ft b + \e ab C, bc Q c (78) 
D. Constraint 



e ab u c R, 



abc • 

a\-\b fn a i\ ( \ ■) r \ 1 i - /-ac\^b 



S a n a +£a b 5 a Z b = (2A - 0) fi - 3£S + e ab ( ac Y l b c + U (79) 



N bc R- a 



abc ■ 

\5a<t> ~ s ab 6 b C ~ S b Cab = -2£e Qb a b - Q (il a + £ Qb S b - 2s ab a b ) - (\6 - |S) (S a - e ab n b ) 

- (s b - e bc n c ) s ab - in a - £ a (80) 

From (C 3 ) ab n b and (Ci) a , or n a u c R abc 

<5 a £ - §<5 Q + 2e ab 5 b n + 25 b S ab = -0 (S a - e ab fl b ) - 2£ (fl a - 3e ab S b ) - 4fie ab ^l b 

+2C ab S b + 2e ab C bc n c + S ab a b - 2e ab H b - Q a (81) 



Finally, we note that the equation formed from {Cj,) ab n a n b is equivalent to Eqs (|103[) and (|T9|) . 

It is worth noting that one of Eqs. I|54[) . ([FT]) is redundant since f|54[) — h (|55|) — ([57)) . Also, note that there are 

no evolution equations for A, A a ,a a , and there is no propagation equation for a a ; these all must be determined by 
specifying a choice of frame. 
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E. Maxwell's Equations 

For completeness we also give the decomposition of Maxwell's equations, previously given in [io| . We decompose 
the electric and magnetic field vectors as 

E a = gn a + g a , (82) 
B a = @n a + M a , (83) 

while the 3-current may be written as 

f = Jn a + J a . (84) 

Maxwell's equations then become: 

g + S a g a = -cf>g + g a a a + 2VL.53 + 2fl a ,^ a + fi p c , (85) 

,% + 5 a ,% a = -4>S$ + ,<% a a a -2flg -2fl a g a , (86) 

g - e ab 5 a @ b = 2^ + g a a a - (§0 - £) g + ^ a g a + e ab (A a ^ b + VL a g b ) - ^ J, (87) 

,% + e ab 8 a g b = -2^ + ^ Q a Q -(§<9-E)^ + E a ^ a -e ab 

Sa + s ab (@ b - 5 b d) = i® a - {\4> + A) e ab @ b - (§0 + iS) g a - Qe ab g b 

+g (-a a + E a + e ab n b ) + ,%e ab (A b - a b ) + ^ ab g b - e ab Q bc S§ c - M o Ja, (89) 

- £ab (g b - 8 b g) = -& a + {\4> + A) e ab g b - (§0 + iS) ^ a - rte ab M b 

(-ota + S a + e ab n b ) - ge ab (A b - a b ) + E ab ^ b + e ab ( bc g c . (90) 

Here, MKS units are used (po)i an d Pe is the charge density. The first two equations arise from the constraint ME, 
while the rest are the evolution ME. In flat space in the absence of currents and charges the rhs of these equations 
vanish (for a static 'natural' choice of frame). Thus, gravity modifies ME in the form of generalised currents. Note 
how the rotation terms f, il and il a flip the parities of the EM fields. 



IV. PERTURBATIONS OF SPHERICALLY SYMMETRIC AND LRS SPACETIMES 

The utility of the approach presented here is that for LRS spactimes, for which all quantities are rotationally 
symmetric about a preferred spatial direction (i.e., they admit a one-dimensional isotropy group), all the non-zero 
1+1+2 variables are scalars. This direction may be specified, for example, by a non-degenerate eigenvector of the 
electric Weyl tensor, or by the vorticity vector. A full discussion of LRS spacetimes in the covariant approach is given 
in [3]; see their Table 1 for a summary of the different cases which can occur, in a notation similar to that presented 
here. 

The fact that background quantities are scalars in LRS spacetimes means that under linear perturbations, all vector 
and tensor quantities are automatically gauge invariant, by the Stewart- Walker Lemma [171 ]. We shall now give an 
overview of how to set up the perturbation equations. 

In the background, which we shall take as a general LRS spacetime, all vector and tensor equations are automatically 
zero, resulting in the set 

= (§0-E) (A- \<t>) +2£n + Q, (91) 
{ = (§E-i0)£+(.4-§0)n+iW, (92) 
ft = +.A£ + (E - §0) , (93) 
A- 6 = -(^ + 0)^+i<9 2 + f£ 2 -2ft 2 + i(/i + 3p)-A, (94) 
£ - f A = | (2A - <j>) A - (§0 + §£) £ - fft 2 - £ + in, (95) 
fi + Q = -Q(ji + p)- (ct> + 2A)Q- |En, (96) 
Q+p + tl = -(%<!> + A) II - (±6 + E)Q-(fj, + p) A, (97) 

£ + in + iQ = + (3 S _0)£_!(!0 + ! S) n 

+i (§0 -2A)Q + 3£H - | (/i + p) E, 
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n 



■V - 

n = 
in = 
w = 

= 



-3££ + (0 + §E) W + HQ + fen, 

-i</> 2 + 2£ 2 + (§0 + s) (fe - s) - 1 (/x + A) - in - e , 

-0£+(|0 + £)fi, 

+ (.a - n, 



in) 



- (3f + n + p - in) n - Qt, 

{2A -4>)Cl- 3£E + W. 



(99) 
(100) 
(101) 
(102) 
(103) 
(104) 
(105) 
(106) 



X 



(107) 



These equations were first presented in this form in [12| for LRS Class II models (which satisfy £ = £1 = 0=>7^ = 0, 
and were shown to be consistent with the commutation relation (|34p . 
It is perhaps easier to think of these in matrix form. Let 

( 4>\ 
e 

£ 
A 

n 

w 
/I 

p 
Q 

vny 

be the column matrix of all non-zero scalar quantities. Depending on the LRS model in question X will not be 
this big. For example, for the Schwarzschild solution we have just X = (cf>,A,£) T . Then, in general, this system of 
equations may be cast in the form 

aX + f3X = TX + X T AX (108) 

where the constant matrices a, /3, T, A may be read off from the above equations. 
We can now set up the perturbative proceedure schematically as follows: 

1. Find a complete set of gauge-invariant perturbation variables. This may be achieved by defining 

*a = S a X ; (109) 

i.e., by taking angular derivatives of the background variables we find a new set of gauge-invariant variables. 
The remaining GI variables are all the 1+1+2 vectors and tensors: Xa — {£ a > a a, ■ ■ -),Xab — {Cab, £ab, T~tab, ■ ■ ■), 
which obey linearised versions of the above 1+1+2 equations. Under perturbations Eq. (|108[) becomes 



aX + px = TX + X T AX + A5 a Xa + Be ab 8 a X b 



(110) 



where the matrices A, B have constant coefficients. Evolution and propagation equations for the new GI vari 
ables \P a may be found by taking the angular derivative of Eq. (|110[) , and using the commutation relations 
and (j36|) , giving: 



a* a + /3* a = [r- 



- \4>) a] * Q - (flat + e/3) erf 
*IM + A5 a S b X b + Be bc S a S b X c 



(111) 



These equations replace the corresponding system (|110[) in the 1+1+2 equations. 



2. Harmonic analysis Two parities of harmonics may be introduced, generalising the axial and polar modes for 
spherical symmetry. These were first defined in 0, [T^ , and are discussed in Appendix [Bj These are analogous 
to the scalar-vector-tensor decomposition in FLRW models. After this, all variables become scalars, which are 
functions of two affine parameters associated with u a and n a . 



11 



3. Master Variables At this stage the governing system of equations is linear in the perturbation variables 3? and 
<&, which are the column vectors containing all the even and odd harmonically decomposed variables, and splits 
into two parities. We then have two linear systems of equations looking like 

7* + A* = H$ (112) 

where S is a matrix with coefficients depending only on the background parameters (as well as the harmonic 
index k), and 7, A are constant matrices. The true degrees of freedom of this system will be governed by a 
reduced set of frame independent master variables, which will obey a closed set of wave equations. Finding 
these can be tricky. All other variables are related to the master variables by quadrature, plus frame degrees of 
freedom. See H for the full details in the Schwarzschild case. 

These are the key steps required given a particular LRS model is chosen. Steps 1 and 2 are algorithmic; step 3 can 
be very difficult. 



V. SUMMARY 



We have presented a new semi-tetrad approach to analysing Einstein's field equations. By introducing a single 
space-like vector into the 1+3 approach we decomposed the 1+3 equations into a system of evolution, propagation 
and constraint equations. These were supplemented by a 1+1+2 decomposition of the Ricci equations for the spatial 
vector. Although presented in restricted form elsewhere, the full system was presented here for the first time. 

A key feature of the approach is that under a complete decomposition all objects are covariantly defined scalars, 
2-vectors in the sheet and transverse-traceless 2-tensors, also in the sheet. In an LRS spacetime, provided the 
spatial vector is chosen appropriately, all the vectors and tensors vanish, leaving just scalars. Under perturbations 
all indexed objects are first-order ensuring that there are no tensorial products; this ensures that we can introduce 
natural harmonic functions on the background which remove all tensorial properties of the equations. Finally, we 
are left with a system of gauge-invariant and covariant first-order PDEs to manipulate. The solution of this system 
provides the solution of the perturbation problem. 
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APPENDIX A: USEFUL RELATIONS FOR DECOMPOSING EQUATIONS 

Given any 1+3 vectors and tensors, we may decompose them as 

x a = Xn a +X a , 
if = Yn a + Y a , 
i>ab = ip(ab) = * {n a n b - \N ab ) + 2^ (a ri b) + * af) , 

\N nh ) + 2$ (a n b) + $ afc . 



4>ab = 4>(ab) = $ (n a rib 2 
Then we have the following expansions from 1+3 quantities 



1+1+2 variables: 



= (e bc X b Y c ) n a + e ah (YX b - XY b ) 



= X 2 +X a X a 

VabcX"y° 

i> ab x b = (X^ + X b ^ b )n a -^X a + X^ a + ^ ab X b , 
Vcd(aX c A) d = e cd X c ^ d (n a n b - \N ab ) + \(X^ C - |*A C ) e c(a + e cd X c ^ d {a 



i (2XY - X C Y C ) (n a n b - \N ab ) + [XY (a + YX (a ] n b) + X {a Y t 



n b ) 



(Al) 
(A2) 
(A3) 
(A4) 



Ipabi) 



a b 



■q, 2 + 2* n * Q 



(A5) 
(A6) 
(A7) 
(A8) 

A c e c{a * b} ,(A9) 
(A10) 
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For 1+3 derivatives we find: 

(X -X b a b ^j n a + Xa a + X a , 



X (a ) 
1p(ab) 

D Q x Q 

VabcV b X C 



- 2* c a c ) n a n b - \^N ab + 
X + Xcb-X a a a + 5 a X a , 
(2XC + e bc S b X c ) i 



3*a ( „ + 2* (a - 2a c * c(a 



n b) 



2* 



(a«b) 



■ (,X a + e a b 



V( a X b) 



3 

+ 



-Xa b + 6 b X - X b - \4>X b - ( bc X, 
2X-cj)X- 2X c a c - <5 c X c j (n a n b - \N ab ) 

Xa {a + 5 {a X + X {& - \4>X {a + X c (£e c(a - ( c(a ) n b) + X( ab + 6 {a X b} , 
(if + - 2* 6 a 6 + <5 b * b - * 6c C 6c ) n a + * a + §0* a + |*a a - ±<5 a * 

-*abCL b + [-CSab + Cab] + ^ ab , 

(3£* + e cd 5 c ^ d - s cd ^ de C e ) (n a n b - \N ab ) 

+ { [-f^* + + + l^c + ^^cdj £c{a + 5 ^ (a + £ «i ^ d(c{a + I nfc) 

-£c{J c ^a} + 2e c{a a c * fc} + £ c{a * c b} + \<t>e c{a ^ c b} - ^e c{a ( c b} + ^ ab + e c{a * b}d C 
Given any relation in 1+3 notation, these relations may be substituted directly to aid decomposition 



cd 



(All) 
(A12) 

(A13) 
(A14) 
(A15) 
(A16) 

(A17) 
(A18) 

(A19) 



APPENDIX B: HARMONIC FUNCTIONS 

We introduce dimensionless harmonic functions Q, defined on any LRS background, as eigcnfunctions of the 2- 
dimensional Laplace-Beltrami operator: 

5 2 Q = -^Q, Q = = Q (0<fc 2 ). (Bl) 
The function r is, up to an irrelevant constant, covariantly defined by 

H'"^' «- r -°- (B2) 

While we haven't chosen a specific basis for Q, we can now expand any first order scalar tp in terms of these functions 
schematically as 

V^E^Q^^sQ, (B3) 

k 

where the sum (or integral) over k is implicit in the last equality. The S subscript reminds us that ip is a scalar, and 
that a harmonic expansion has been made. 

We also need to expand vectors in harmonics. We therefore define the even (electric) parity vector harmonics as 

Q a k) =rS a Q^ => Q a = = Q a , 5 2 Q a ={l-k 2 )r- 2 Q a ; (B4) 

where the (k) superscript is implicit, and we define odd (magnetic) parity vector harmonics as 

Qi k) = re ab 6 b QW Q_=o = Q_, s 2 Q a = (l - k 2 ) r - 2 Q a . (B5) 

Note that Q a = e ab Q b <*=> Q a = —e ab Q b , so that e ab is a parity operator. The crucial difference between these two 
types of vector harmonics is that Q a is solenoidal, so 

S a Q a = 0, (B6) 
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while 

5 a Q a = ~k 2 r- x Q. (B7) 

Note also that 

e ab S a Q b = 0, and e ab S a Q b = +k 2 r - 1 Q. (B8) 

The harmonics are orthogonal: Q a Q a — (for each k), which implies that any first-order vector ip a can now be 
written 

= ^ k) Q ( a } + ^ k) Q { a k) = ^Qa + Ma- (B9) 

k 

Again, we implicitly assume a sum over k in the last equality, and the V subscript reminds us that ip a is a vector 
expanded in harmonics. 

Similarly we define even and odd tensor spherical harmonics as 

Qab = r 2 5 [a S b} Q, =*> Q ab = Q = Q ab , (BIO) 
Qab = r 2 e c{a 6 c 5 b} Q, => Q ab = = Q ab , (Bll) 

which are orthogonal: Q ab Q ab = 0, and are parity inversions of one another: Q ab = —e c { a Q b ^ Q ab = £ c {aQ b f- 
Any first-order tensor may be expanded 



*ab = J2 ^Qab + *T ^ab = ^Qab + ^jQab- (B12) 

k 
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